Dielectric Function of Diluted Magnetic Semiconductors in the Infrared Regime 
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We present a study of the dielectric function of metallic (III,Mn)V diluted magnetic semicon- 
ductors in the infrared regime. Our theoretical approach is based on the kinetic exchange model 
for carrier induced (III,Mn)V ferromagnetism. The dielectric function is calculated within the ran- 
dom phase approximation and, within this metallic regime, we treat disorder effects perturbatively 
and thermal effects within the mean field approximation. We also discuss the implications of this 
calculations on carrier concentration measurements from the optical f-sum rule and the analysis of 
plasmon-phonon coupled modes in Raman spectra. 
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I. INTRODUCTION 

Ever since post-growth annealing procedures in 
(III,Mn)V diluted magnetic semiconductors (DMS) 1 
demonstrated the ability to increase the ferromagnetic 
transition temperature, T c , by almost two fold, the study 
of these promising and interesting materials has been 
vigourous both experimental and theoretical. The ob- 
served increase in T c is due to the increase of the free hole- 
carrier concentration originally induced by the Mn sub- 
stitutional impurities but partially compensated by other 
impurities such as As-antisites and Mn-interstitials which 
arise from the needed non-equilibrium growth of these 
materials. 2-6 Because of the sensitivity of their mag- 
netic state to growth conditions, doping, and external 
fields, and the strong valence-band spin-orbit coupling 
the transport and optical properties of these heavily- 
doped semiconductors are richer than those of conven- 
tional itinerant electron ferromagnets. 

We present a theory of the infrared dielectric func- 
tion of (III,Mn)V ferromagnets based on the kinetic ex- 
change model for carrier induced (III,Mn)V ferromag- 
netism which assumes a sallow acceptor picture. This 
justifies the treatment of the electronic structure of the 
free carriers by the unperturbed host semiconductor band 
structure. The Mn 3d 5 electrons are strongly local- 
ized with a fully polarized 5 = 5/2 moment which in- 
teracts with the free carriers via Coulomb and short- 
range exchange interactions, 7-9 as has been demon- 
strated by electron paramagnetic resonance (EPR) and 
optical measurements. 2,10 As described in the next sec- 
tion, the dielectric function is calculated within the ran- 
dom phase approximation and, in this metallic regime, 
we treat disorder effects perturbatively and thermal ef- 
fects within the mean field approximation. 

In particular, we demonstrate that infrared dielec- 
tric function measurements in metallic samples can be 
used to obtain information about the carrier concentra- 
tion, to simulate the plasmon-phonon overdamped cou- 
pled modes, and to explain naturally the observed optical 
absorption measurements in the infrared regime arising 



from inter-valence-band transitions. Some of our con- 
siderations are based on standard linear-response theory 
expressions for weakly disordered metals, in which dis- 
order is included through finite quasiparticle lifetimes 
and localization effects (demonstrated in this model 
through finite-size exact diagonalization studies 11 ) are 
taken into account phenomenologically by weighting dif- 
ferently intra-band and inter-band quasiparticle life- 
times. We estimate the magnitude of these lifetimes 
from extensive theoretical studies in dc-transport stud- 
ies which describe successfully the metallic DMS high-Tc 

1 2 1 3 

regime. ' 

Although the predictions discussed below are intended 
to be more reliable for the most metallic systems 14 , they 
appear to explain much currently available infrared opti- 
cal data, many of which has been obtained in studies of 
systems with relatively low dc conductivities. 

The paper is organized as follows. In Sec. II we briefly 
describe the model Hamiltonian and our theoretical ap- 
proach and approximations. In Sees. Ill we analyze the 
results of the dielectric function calculations, sum rules, 
and plasmon-phonon coupled modes. We then present a 
summary in Sec. IV 



II. MODEL HAMILTONIAN 

Focusing on the infrared regime, in our calculations we 
describe the free carriers electronic structure of the DMS 
using the six band Kohn-Luttinger Hamiltonian. The di- 
electric function is calculated in the Random Phase Ap- 
proximation (RPA) formalism and finite disorder effects 
are taken into account perturbatively by introducing a 
finite quasiparticle lifetime. In addition, thermal fluctu- 
ations are ignored and finite temperature effects are only 
considered at the mean field level. Hence, the system is 
described by the following Hamiltonian, 

H = H\ 10 \ cs + H Mn 2+ _hoic S + J ^ Sj-Sj S(rj — Rj) , (I) 

I,i 

Where Hhoies is the part of the Hamiltonian which 
describes the itinerant holes within the k • p theory, 
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iJ Mn 2+ _hoios describes the screened Coulomb interactions 
between the holes and localized Mn ions, and the last 
term is the antifcrromagnctic exchange interaction with 
strength J between the spin of the Mn 2+ ions located 
at the random positions R/ and the spin Sj of the hole 
carriers. 

The sources of disorder known to be relevant in these 
materials include the positional randomness, R/ of the 
substitutional Mn ions with charge — e, random place- 
ment of interstitial Mn ions, acting as non-ferromagnetic 
double donors, 15 and non-magnetic As anti-sites acting 
also as having charge +2e. These Coulomb interactions 
are included in -ffMn 2 +-hoics- Previous estimates of the 
valence band quasiparticles lifetimes using Fermi's golden 
rule induced by this term are of the order of 100-250 
meV. 12 We treat H Mn 2+ _ holes within the Born approxi- 
mation in the Kubo linear response formalism as in Ref. 
12,13,16. Hence, in what follows, we drop i?Mn 2 +-hoics 
and introduce its effects through a finite quasiparticle 
lifetime. 

In order to obtain the temperature dependence of the 
electronic and magnetic properties of the system, we min- 
imize the free energy per unit volume, 



3~ 3~i oris -\- J~hole 



(2) 



where Ti ons is the contribution of the ion spins to the free 
energy and in the mean field description has the form, 



-TN Mn 2+ In 



sinh (hS/T) 
(hS/2T) 



(3) 



where h — Jp£/2, p and N Mn 2+ are the free carrier den- 
sity and Mn density respectively, and £ is the spin po- 
larization of the carriers. Tholes is the free energy of the 
holes, which is obtained in the virtual crystal approxima- 
tion (VCA) using a Luttinger k • p model for describing 
the carriers. In the VCA the average density of states for 
the real system is replaced by that of the average Hamil- 
tonian. This approach implies a translational invariant 
system with an effective exchange field acting on the car- 
rier spins h ex = mN Mn 2+ SJ , where m is the polarization 
of the Mn spins. In the Luttinger k • p model the wave 
function of the holes in the state (n, k), where n is the 
subband index and k is the wave vector, is expressed as 



^,k(r)=e ik - r ^a^r|J,m J ) 



J.nij 



(4) 



where \J,mj) are the six Y± v valence band wave func- 
tions. The coefficients a^k'' an( ^ tne corresponding 
eigenvalues, e nj k, depend on the spin polarization of 
the Mn, and are obtained, from the Luttinger k • p 
Hamiltonian. 8,17 

By minimizing Eq.(2) we obtain the T-dependence of 
the spin polarization of the carriers, £(T) and of the Mn, 
m(T). In the VCA the Curie temperature has the expres- 
sion T c — N Mn 2+ S 2 J 2 / 3\p where \p is the zero wavevec- 
tor paramagnetic susceptibility of the hole gas. Through 



the dependence of the Hamiltonian on the spin polariza- 
tion, we obtain the temperature dependence of the eigen- 
vectors, eigenvalues and chemical potential. Through- 
out this paper we consider always the typical optimal 
Mn concentration of x = 0.05 and an exchange coupling 
J = 60eVnm 3 . 7,18 Although the properties of the sys- 
tem depend on the orientation of the magnetization, this 
dependence is much smaller than the T-dependence and 
since we are interested in the variation of the electronic 
properties with T we fix the Mn spin polarization in the 
^-direction. Experimentally this can be easily achieved 
by a small coercive field of the order of 100 Oz. 

From the eigenvalues and eigenvectors, we calculate 
the dielectric function that in the RPA has the form, 



e(q,w) 47re 2 

= 1 + J *( q ' W ) i 



(5) 



being x(q, uo) the susceptibility, 

"F(£j,k+ q ) - n F (£j,k) 

flUJ + Ejfc - £i,k+q 



x(q,w) = 



/„-(k,k + q) , 



(6) 



with 



/ ilj (k, k ')= £Krr«#? j i ■ (") 

where rip denotes the Fermi-Dirac distribution and 
£oc the dielectric function of the host semiconductor 
(£00=10.90 for GaAs). In the limit q = |q| — > 0, the over- 
lap function becomes /i,j(k, k+q) — > 6ij+Aij(k)q 2 , and 
Eq. (5) can be rewritten as: 



e(w) u 2 Aire 2 

= 1 1 H Xinter(u) 

£00 ^ ^no 



(8) 



The second term in the right hand side of Eq. (8) is the 
intra-band contribution and constitutes the Drude weight 
present in dc measurements and it is closely related to the 
plasmon frequency in the metallic regime, 
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Aire 2 y-^ de lM 
e °° i,k 



dk 



(9) 



The third term is the inter-band contribution and has 
the form, 

XmterM = > , J —A^(k) . (10) 

.... nk} + Ei'.k — £iM 

We also consider the coupling of the holes with longitudi- 
nal optical (LO) phonons. Within the RPA approxima- 
tion this can be done by including a phonon term in the 
susceptibility such that the final dielectric function is: 



e(w) u 2 4-rre 2 
= + 

£oo e^o 



Xinter 



,.2 , ,2 
0J TO - u LO 



u) 2 - uj\ + iuJj pho 



(11) 
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where Hujto — 33.25meV, hiu^o — 36.23meV. In the 
case of the coupled plasmon-phonon systems, the damp- 
ing of the optical phonons is small compared to that of 
plasmons, 19 and it is estimated to be ^phonon /uj TO « 0.01 
for GaAs. We account for disorder by including the life- 
time broadening of quasiparticle spectral functions in a 
phenomenological way. We use different lifetimes for the 
inter and intra-band contributions to the dielectric func- 
tion. In general we expect the intraband lifetime 1/ju 
to be much smaller than the interband lifetime l/7y as 
the intraband scattering occurs at zero energy whereas 
the interband scattering occurs at finite energy This 
has been demonstrated by finite size exact diagonaliza- 
tion studies where a much larger intra-band disorder re- 
duction was observed than the interband contribution at 
finite frequencies. 11 In the previous formalism we have 
considered only the diagonal component of the dielectric 
function parallel to the carrier spin polarization. 



III. RESULTS 
A. Dielectric function 

We have computed the dielectric function for differ- 
ent electron densities, lifetimes and temperatures. We 
illustrate the low-temperature behavior of the dielectric 
function and the ac conductivity, <x(c<j), in Fig. 1 where 
4:Tt<j(uj) — wlm[e(tj)] and — Inr[^y] are plotted for differ- 
ent degree of disorder 7 = 7,, = jij . The conductivity, 
which is finite in the dc limit, has different high-frequency 
contributions which can be explained in terms of inter- 
band transitions: in the clean limit, at this carrier den- 
sity, there are four strong resonances at ui ~ 100, 175, 
275 and 390meV, due to transitions between the spin 
polarized six k • p hole bands. Upon increasing 7, the 
Drude peak contribution (intraband) broadens and over- 
laps with the high-frequency interband contribution. Our 
results agree with those obtained previously by Sinova et 
al in Ref. 16 and are consistent with the experiments of 
Singley et al in Ref. 20 and 21. Nonetheless, the dc 
conductivity in the experiments is always much smaller 
than the ac conductivity around li ~ 2000cm -1 . This 
is an indication that intraband disorder is larger than 
interband disorder in the experiments. This is also in 
agreement with finite size exact diagonalization studies 
where a much larger intra-band spectral weight reduc- 
tion was observed due to disorder induced localization. 11 
We illustrate this in Fig. 2a where we plot 47rcr(w) vs 
lu for fixed interband disorder 7^ =25 meV and different 
intraband disorder 7^. For 7^ >> 7^ the dc limit is 
smaller than the ac conductivity for frequencies up to 
u> > 2000 cm -1 (for 7^=250 meV) which corroborates 
our previous argument. 

In Fig. 3 we plot 47rlm[tr(a;)] and — lm[l/e(w)] for dif- 
ferent holes densities in the case of zero temperature, 
7ii=100 meV, and 7^=25 meV. As obtained by Sinova 
et al in Ref. 16, we find that , as the hole density increases 
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FIG. 1: (Color online.) Frequency dependence of the dielec- 
tric function, a) 4tv<j(u) — wlm[e(a;)]. b) — Im[^4jj]. c) 
Same as in b) showing a blow-up of the coupled plasmon- 
phonon mode region. Parameters of the calculation: T = OK, 
p — 0.5nm -3 and different interband and intraband scattering 
times: 7,; = 7^=2.5, 25, 50, 100 and 250meV. 



a peak at energy near 280 meV emerges and the Drude 
peak appears more clearly. 



The main features of — Im[^ijy](Fig. 3b) are a strong 
single plasmon-likc peak (marked in the figure with a cir- 
cle for p=0.2 nm~ 3 ) with much weaker shoulders (marked 
with squares) that correspond to interband transitions. 
The plasmon-like character of the strong peak is evident 
when we plot the position of the peaks in — Im[^jy] as 
function of the hole density (Fig. 4). The energy of the 
strong peaks (circles) can be fitted perfectly to the ex- 



pression 



with to p /=0.25 m e . The energies of 
~230 meV and ~375 meV (squares) 



47re 2 p 

'-<x>m pl ' 

the other peaks at 
depend very little on the hole density, and we assign them 
to the light-hole band to heavy-hole band transition and 
to the light-hole band to the split-off spin-orbit band 
transition respectively. Of course these transitions are 
coupled with the plasmon-like excitation. Note that the 
peak positions are also rather independent of the carrier 
spin polarization, as the concentration of Mn impurities 
is fixed at x=5% for all the hole densities, and the carrier 
polarization changes from 0.81 for p=0.1nm~ 3 to 0.57 for 
p=0.6nm -3 
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FIG. 2: (Color online.) Frequency dependence of the dielec- 
tric function, a) ulm[e(u)] = 4ira(u>). b) — Im[j7^y]. c) 
Same as in b) showing a blow-up of the coupled plasmon- 
phonon mode region. Parameters of the calculation: T = OK, 
p — 0.5nm -3 , fixed jij=25 meV and different ju: 2.5, 25, 50, 
100 and 250meV. 
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FIG. 4: Energy positions of the peaks appearing in Im[ Jjjy] 
as a function of the hole density p. Circles corresponds to the 
stronger peak in the spectra whereas square dot correspond to 
the weaker peaks. The dashed line is the energy of a plasmon 
with an effective mass 0.25 times the free electron mass. 



B. Sum rules 
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We can also define an optical effective mass from the 
sum rules of the model as 
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FIG. 3: (Color online) Frequency dependence of the dielectric 
function, a) wlm[e(w)] = 4iva(uj). b) — Imfjr^j]. c) Same 
as in b) showing a blow-up of the coupled plasmon-phonon 
mode region. Parameters of the calculation: T — OK, 7^ =25 
meV, 7;i = 100meV, and hole densities varying from O.lnm -3 
to 0.5nm -3 . The symbols mark the position of the plasmon- 
like peak (circle) and of the peaks corresponding to interband 
transitions (squares) for p=0.2nm -3 (see text). 
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Both sum rules, when integrated up to infinity give the 
same optical mass. However, because the calculation can 
only be done in a finite frequency range, in practice the 
two sum rules give different optical masses that depend 
on the frequency cut-off in the integration, w cuto ff. In 
Fig. 5, we plot the quantity p/m op t as obtained from 
the two different sum rules as a function of cj cu toff for 
p=0.5nm -3 . Note that it is necessary to sum up to rather 
large values of the frequency in order to get the same ef- 
fective mass from the two different sum rules. Also, we 
see that in the six-band k • p Hamiltonian it is neces- 
sary to integrate up to very high energies for obtaining 
an accurate value of p/m opt . That imposes a serious lim- 
itation for the experimental accurate estimation of this 
quantity to a few percent but can still be used in com- 
bination with other indirect transport measurements to 
establish consistency. 

Recently, Singley et aP 1 have studied experimentally 
the dependence of p/m opt on the cutoff frequency and 
temperature by integrating the measured infrared con- 
ductivity up to 1.5 eV. In order to compare with these 
experiments, we study p/m op t as obtained from the con- 
ductivity for different hole densities. The value of the 
optical mass so obtained depends on the density of holes 
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FIG. 5: Cut-off frequency dependence of the two sum rules. 
The results correspond to T = OK, p — 0.5nm~ 3 ,ju — 
lOOmeV and 7y = 25meV. 



and, for a finite frequency cutoff, on the lifetimes used in 
the calculations. This is illustrated in Fig. 6 where we plot 
the conductivity sum rule as function of this frequency 
cut-off for different temperatures, densities and quasipar- 
ticle lifetimes. ^From the figure we see that in the most 
metallic samples, lowest "fa, the sum rule will be most 
accurate in determining the carrier concentration. How- 
ever, it is clear that for moderate hole densities the value 
of p/m op t obtained from the sum rule is far from the con- 
verged value. In any case the values of the optical mass 
obtained from the sum rule are in the range O.I7-0.30m e 
depending on the density, frequency cutoff and temper- 
ature. In Fig. 6 we also compare the value of p/m op t at 
zero temperature (ferromagnetic phase) with its value at 
high temperature (paramagnetic phase). At low temper- 
atures there is an enhancement of the low energy spectral 
weight. This can be attributed to low energy interband 
transitions in the spin polarized system which originate 
from the spin-orbit coupling. These transitions are ab- 
sent in the paramagnetic phase because there is not a 
Zeeman-like spin splitting of the bands. In the absence of 
spin-orbit coupling, these excitations correspond to tran- 
sitions between bands with well defined spin antiparallel 
polarization which, therefore, do not contribute to the 
charge-charge response function, and are absent in the 
Im[^ijj] spectrum. 



C. Plasmon-phonon coupled modes 

For the range of carrier densities considered, a 
plasmon-phonon coupled mode occurs at lo w u>to- This 
plasmon-phonon coupled mode (Fig. lc) is very sensitive 
to disorder: it changes from a peak-like shape to a Fano- 
like shape upon increasing 7^. Interestingly, only intra- 
band disorder modifies significantly the coupled mode: 
-Imf^rr] around u> « loto changes little when 7^ in- 
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(Color online) Cut-off frequency dependence of the 
: ^ dwReKw)]. a) 7 « = 7lj = 
OK (black solid line) and T = 200K (red 
dashed line), for p — O.lnm -3 (bottom curves) and p = 
0.5nm -3 (top curves) b) Same as in a) but ju = 250meV 
and 7y = 25meV. 



creases. This can be seen by comparing the curves in 
Fig. lc (where 7^ varies from 2.5meV to 500meV) with 
the corresponding curves (i.e., with the same 7^) in 
Fig. 2c where 7y=25meV is fixed. Due to the energy 
difference between ujlo an d the interband transition en- 
ergies, the phonon modes almost do not couple with the 
interband transitions. The plasmon-like part of the cou- 
pled mode is superimposed to the plasmon peak for small 
7«. On the other hand, for ja > 250 meV the plas- 
mon peak is strongly overdamped (this corresponds to 
a phonon peak at u = ojto superimposed to an almost 
featureless a(u>), see Fig. la) and only the phonon-like 
coupled mode remains. This is likely the situation cor- 
responding to the Raman experiments in Ref. 22 and 
23. Indeed, the authors of Ref. 23 model their data by 
using a simplified Drude model where interband transi- 
tions are neglected, and obtain the value of the carrier 
density indirectly by a line-shape analysis of the coupled 
plasmon LO phonon in their Raman spectra. By using 
an averaged hole mass, they conclude that the same line- 
shape is obtained if one keeps the ratio ^- fixed, which 
allows them to estimate the carrier density for different 
values of the plasmon damping. Note, however, that in 
a simplified Drude model s(u>) oc — P tT . which, together 
with our results, suggests that it should not be possible 
to know p, m opt and 7^ independently. Thus, we be- 
lieve that further experimental work is needed in order 
to clarify this important point. 
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IV. SUMMARY 

In closing, we have presented a theory of the infrared 
dielectric function of (III,Mn)V ferromagnets based on 
the kinetic exchange model for carrier induced (III,Mn)V 
ferromagnetism. Our results demonstrate that the in- 
frared dielectric function in metallic samples can be used 
to obtain information about the carrier concentration, 
to simulate the plasmon-phonon overdamped coupled 
modes, and to explain naturally the observed optical 
absorption measurements in the infrared regime aris- 
ing from inter-valence-band transitions. The predicted 
p/^opt ratios could be tested in the metallic regime 
within the newly grown samples which exhibit large con- 
ductivities. The utilization of the f-sum rule to estimate 
the carrier concentration has been shown to be most 
accurate within the highly metallic samples and should 
be tested experimentally by careful comparison of these 



measurements with other means of measuring the car- 
rier concentration. We also have shown that the analysis 
of the Raman plasmon-phonon coupled modes within a 
simple one band Drude model as in Rcf. 23 arc not very 
accurate since estimates of m opt and ju cannot be done 
independently 
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